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ABSTRACT 

The galaxy M51 was observed using the Mimir instrument on the Perkins telescope to constrain 
the resolved H-band (1.6 /im) polarization across the galaxy. These observations place an upper limit 
of Ph < 0.05% on the i^-band polarization across the face of M51, at 0.6 arcsecond pixel sampling. 
Even with smoothing to coarser angular resolutions, to reduce polarization uncertainty, the i^-band 
polarization remains undetected. The polarization upper limit at i^-band, when combined with pre- 
vious resolved optical polarimetry, rules out a Serkowski-like polarization dependence on wavelength. 
Other polarization mechanisms cannot account for the observed polarization ratio (Ph/PVRI < 0.05) 
across the face of M51. 

Subject headings: Galaxies: individual: (M51) — Galaxies: ISM — Galaxies: magnetic fields — 
Infrared: galaxies — Polarization — Radiation mechanisms: general 



1. INTRODUCTION 

A comparative study of magnetic fields probed at op- 
tical, near-infrared (NIR), and radio wavelengths for a 
face-on galaxy can elucidate the nature of galaxy-scale 
magnetic fields and provide context for studies of mag- 
netic fields in the Milky Way. To this end, the galaxy 
M51 was chosen for a pilot study of resolved pol arimetry 
with the Mimir instrument (j Clemens et al .1120071 ). M51 is 
an essentially face-on, grand-design spiral galaxy, though 
not a perfect Mil ky Way analog because it is interacting 
with NGC 5195 (jToomre fc Toomrd[l972h . Its high ele- 
vation (declination +47°) made it suitable for multiple, 
long observations and its large angular size (~ 11 x 7 
arcmin) was well matched to the Mimir 10 x 10 arcmin 
field-of-view. 

How does the magnetic pitch angle across the face 
of M51 vary with position at NIR wavelengths, and 
how does this compar e to pitch angle va riations seen 
at other wavelengths (Patrik eev et al.l l20Q6h ? How do 
these pitc h angles compare to results from the Milky 
Way f e.g.,|Han_fc Qiao 1991; lHeiledll99?l iPshirkov et all 
l20ia iPavel et al.l I2012D ? Since optical/NIR polariza- 
tion traces magnetic fields in dusty regions while syn- 
chrotron polarization and Faraday rotation trace mag- 
netic fields in the warm interstellar medium, are the 
same large-scale magnetic fields being observed in both 
regions, or are local magnetic fields able to decouple from 
the large-scale field under certain circumstances? Does 
the resolved polarization across the face of M51 follow 
a Serkowski dependence of polar ization strength with 
wavelength ( Serk owski et al.l [1975) as seen for Galactic 
lines-of-sight? 

The magnetic field of M51 has been previously stud- 
ied at optical and radio wavelengths, but never at NIR 
wavelengths. The first observations of th e magn etic field 
in M51 were made bv iMathewson et al.l (|1972f ) at 1415 
MHz (« 21 cm) using synchrotron emission. Later, 
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iSegalovitz et al.l (jl976h observed polarized synchrotron 
emission at 6 and 21 cm, but found no measurable polar- 
ization at 49 cm. The observed polarizations suggested 
that the large-scale magnetic field is oriented parallel to 
the spiral arms of M51. A dditional polarimet ric obser- 
vations by iNeininge 3 (|1992l) at 2.8 cm iHorellou et al.l 
(|1992f ) at 1 8 and 20.5 cm iHeald et al.l (|2009f ) at 18 and 
22 cm, and lFletcher et all ( 2011f ) at 3 and 6 cm support 
this interpretation. However, in all cases, those authors 
warn that Faraday depolarization may affect interpreta- 
tion of the observations. 

To probe the magnetic field of M5 1 without t he ef- 
fects of Faraday depolarization, Sca rrott et al.l ([1987) 
used resolved optical polarimetry (Pvri- 450-1000 nm, 
with peak response at 850 nm) to measure the orien- 
tation of disk magnetic fields (as projected onto the 
plane of the sky) assuming the polarization was caused 
by magnetically- aligned dust grains. Those observations 
also traced an open spiral pattern from the nucleus to a 
radial offset from the core of 4-5 kpc. 

This pilot s tudy intended to extend the work of 
iScarrott et al.l (|1987l ) by measuring resolved polarimetry 
across the disk of M51 in the NIR. However, the fully re- 
duced H-band polarimetric observations of M51 showed 
a lack of polarization. After coadding multiple obser- 
vations, a very low upper limit on the H-band polariza- 
tion was set. A number of polarization mechanisms are 
considered, but none can account for the observed and 
absent polarizations. 

2. OBSERVATIONS AND ANALYSIS 

M51 was observed in 2009 February, March, and May 
with Mimir on the Perkins 1.8m Telescope outside of 
Flagstaff, AZ. Fifty observations were obtained, totaling 
14.5 hours of on-target integration. Imaging polarime- 
try was measured by a combination of a cold, stepping 
half- wave plate (HWP, to modulate the polarization sig- 
nal) and a cold, fixed wire grid. For each pointing, 
the HWP was rotated to 16 unique position angles (the 
first HWP position was measured twice per HWP rota- 
tion, for a total of 17 images). At each HWP angle, a 
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Fig. 1. — (Left Panel) Linear display of the Mimir 9x9 arcmin H-band photometric image of M51. This image is the result of the 
coadding of all observations of M51. (Middle panel) Measured polarization before statistical debiasing (P') at the native 0.6 pixel sampling. 
(Right panel) Measured polarization uncertainty across M51 at the same sampling. 



10.25 sec exposure was taken. This procedure was re- 
peated toward six different sky dither positions, offset 
by 15 arcsec, following a rotated hexagon pattern on the 
sky, for a total of 102 images per pointing. Resolved 
polarimetry with Mimir in this manner has successfully 
be en used for the c omet 73P / Schwassmann- Wachmann 
3 (| Jones et al.ll2008[) and for the reflection nebula IRAS 
05329-0505 (jPavell l2013h . The Mimir data reduction 
pa ckages (developed fo r starlight polarimetry) described 
in lClemens et al.l (2012b.c) were applied to these obser- 
vations, but with the key differences described below. 

2.1. Background Fitting 

To remove sky backgrounds, the Mimir reduction soft- 
ware normally calculates 'super-sky' images for each 
HWP position after flat-fielding, to remove the second- 
order (illumination and other) effects. To do so, for each 
HWP position during a single observation, all images 
taken through that HWP position are coadded without 
astrometric registration to reject stars and the result is 
subtracted from the constituent HWP images. Since the 
sky background may be polarized, it is important to pre- 
serve this polarization information by calculating such 
sky backgrounds independently for each HWP angle. For 
M51, super-sky images created in this manner would end 
up removing much of the object's flux (since the 15 arcsec 
dither size is smaller than M51). Instead, backgrounds 
were modeled and removed using a different method. 

For each pointing, 'Sky only' zones were identified by 
eye in the Mimir images of M51 and used to fit two- 
dimensional, second-order polynomials to the sky back- 
ground across each image. While these zones were iden- 
tical in all 102 images per observation, the fitted sky 
background polynomials were independently fit in each 
image. These zones avoided M51, NGC 5195, and the 
brightest stars in the field. The choice of a second-order 
polynomial was based on an analysis of the sky back- 
ground residuals. These polynomial fits were then sub- 
tracted from the images to remove the sky background 
prior to astrometric registration and co-adding. 

2.2. Resolved Polarimetry 

The analysis of reso lved polarimetry is sim ilar to the 
analysis described in iClemens et al.l (|2012d ) for unre- 
solved starlight polarimetry. After removing the sky 



background, images were grouped by HWP angle and 
coadded to form 16 master HWP images for each obser- 
vation. Since the observations were dithered, this step 
removed image defects (e.g., bad pixels and cosmic rays). 
Because of the 40 modulation of the polarization as the 
HWP rotated, the 16 master HWP angles correspond to 
four independent measurements of each of the instrument 
independent polarimetric position angles (IPPAs): 0, 45, 
90, and 135°. The four images at each IPPA were aver- 
aged to form four master IPPA images per observation. 

From these four master IPPA images, corrected Stokes 
images (Q/Ucorr) were calculated: 
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^CORR = (^45 - ^135 )/ (^45 + ^135) _ ^INST, (2) 

where operations were performed on each pixel in the 
entire image and Q/Ltnst are the previous ly determined 
instrumental polarization contributions (Cle mens et al.l 
I2012b[ ). 

These corrected images were scaled by the polarization 
efficiency (77 = 91.1±0.4%) and rotated by the HWP zero- 
phase offset angle (0, the angle b etween the instrumen t 
coordinate system and equatorial; IClemens et a l. 2012b): 

Q = cos(20)Qcorr/v ~ sin(20)U C ORR/v (3) 

U = cos(20)U C orr/v + sin(20)U C OBR/v- (4) 

This generated a total of 50 unique Stokes U and 50 
Stokes Q images for M51. These were astrometrically 
registered and coadded to produce the final Stokes U 
and Q images for M51. From these corrected, calibrated 
Stokes images, an image of the biased degree of polariza- 
tion (P') was calculated: 



P' = VU 2 + Q 2 



(5) 



Since the degree of polarization suffers from a positive 
statistical bias, P' was deb iased using the prescription of 
IWardle fc Kronbergl (|1974f ). to calculate the true polar- 
ization: 



p/2 



(6) 



A deep photometric image was also calculated by reg- 
istering and coadding all of the individual observations 
and is shown in the left panel of Figure [U 
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3. RESULTS 

The middle panel of Figure Q] shows the biased polar- 
ization percentage at Mimir's native 0.6 arcsecond pixel 
sampling before statistical debiasing and before any at- 
tempt was made to increase polarization signal-to-noise 
ratio (SNR) via smoothing. After applying the correction 
for the positive statistical bias, the polarization unceit 
tainties (Fig. [TJ right panel) are greater than the me^ 
sured polarizations across the entire image. Therefore^ 
there is no measurable H-band polarization across M5% 
or NGC 5195. A polarization position angle image w$§ 
also calculated, but is meaningless when there is no mea- 
surable polarization and so is not shown. 

The uncertainty in the measured degree of polariza- 
tion, at 0.6 arcsecond sampling, is shown in the right 
panel of Figure HJ The polarization uncertainty is almost 
perfectly anticorrelated with the total flux, shown in the 
left panel of Fig. HJ This is expected, since the polariza- 
tion uncertainty depends directly on the uncertainty in 
the flux measured at each HWP angle. If photon noise 
dominates, then the photometric SNR should increase 
with total flux. Additional sources of uncertainty include 
the background subtraction, a conservative 0.1% instru- 
mental polarization calibration uncertainty, and the un- 
certainty in the polarization efficiency (1 part in 225). 

Since the Stokes U and Q parameters are Gaussian- 
distributed quantities, the polarimetric SNR across the 
face of M51 can be increased by smoothing Stokes U and 
Q images to synthesize lower angular resolution images. 
The effects of smoothing on debiased polarization (P) 
and its uncertainty (<rp) across the face of M51 are shown 
in Figure [2] for the median of pixels with emission above 
the mean sky background. Smoothing to 5.8 arcsecond 
resolution (10 pixels) caused the polarization uncertainty 
to decrease and the debiased polarization to increase. 
Despite this, the median debiased degree of polarization 
remained consistent with no detection. Smoothing to 
even coarser resolutions (15-60 arcsec) resulted in larger 
polarizations, but also larger corresponding uncertain- 
ties. Smoothing beyond 60 arcsec caused both the po- 
larization and uncertainty to fall, but the SNR never 
exceeded two. 

Smoothing to such large resolutions will tend to mask 
any bona fide polarization signal because of depolariza- 
tion. If the spiral-like pattern seen by iScarrott et al.l 
(1987) exists at H-band, these last large smoothing ker- 
nels would have smoothed over regions with significant 
rotation in PA and would effectively depolarize the syn- 
thetic pixels. This is critical in the core of M51, where a 
large synthetic pixel would encompass all possible PAs. 
Therefore, the marginal polarization detections at reso- 
lutions larger than 60 arcseconds are likely not reliably 
revealing polarizations across the face of M51. 

4. DISCUSSION 

These results place an H-band polarization la upper 
limit of 0.05% across the face of M51, at an angular 
sampling of 0.6 arcseconds and 0.02% at a resolution 
of 1.8 arcseconds. Low NIR polarization s are expected, 
given the ~ 1% optical polarizations from IScarrott et aLl 
(|1987[ ). If the polarization is caused by dichroic dust 
extinction in the disk of M51, then the wavelength de- 
pendence of polarization should follow a Serkowski law 
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Fig. 2. — Runs of median polarization percentage (P), and its 
uncertainty(crp), versus smoothing resolution, for all M51 pixels 
showing emission brighter than the sky background. Median po- 
larization is shown as the solid curve; uncertainty as the dashed 
curve. The lower x-axis represents the number of pixels (n) com- 
bined into a new (nxn) synthetic pixel. The upper x-axis is the cor- 
responding angular resolution of the new synthetic pixel. Stronger 
smoothing increases the median polarization, but also increases its 
uncertainty, thereby failing to improve the ratio of those quantities. 
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where A max is the wavelength of maximum polarization. 
Assu ming A ma , x = 0.55 /im ( typical for Galactic lines-of- 
sight: ISerkowski et al.l 119751 ) would predict Pn/Pvm — 
0.39, whereas the observations show PH/Pym < 0.05 ± 
0.05. 

To test whether any Serkowski Law could account for 
these observations, the assumption of A max = 0.55 /im 
was relaxed and the NIR and optical polarization data 
were used to fit for the best A max , using the M PFITFUN 
routi ne from the IDL-based MPFIT package (|Markwardtl 
2009). No value of A max could recreate the observed po- 
larization ratio. Polarization by dust dichroism, which is 
obser ved throughout the Milky Way (jMathewson fc Fordl 
1970), cannot explain these observations. 

Other possible polarization sources were considered to 
try to account for the extra optical polarization, includ- 
ing electron (Thompson) scattering of low energy pho- 
tons. However that polarization mechanism is wave- 
length independent (jAntonucci fc Millerl fl985). so any 
observed wavelength dependence would only reflect the 
polarization spectrum of the photon source. The signif- 
icant observed difference between the degree of polar- 
ization at optical and NIR wavelengths rules out elec- 
tron scattering, since it would only enhance the degree 
of polarization, not change Pu/Pvm- Polarized syn- 
chrotron radiation is also unable to account for the ob- 
served polarization ratio, since the degree of polariza- 
tion only depends on the electron power law distribution. 
The synchrotron intensity ratio at optical and H-band 
would have to be at least l%/0.05% = 20. This requires 
either optically thick synchrotron emission or optically 
thin emission with an inverted electron energy distribu- 
tion (both are unlikely). 
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For now, the cause of the observed polarization at op- 
tical and the lack of polarization at NIR wavelengths 
across the face of M51 remains unknown. Resolved po- 
larimetry in multiple passbands may provide additional 
insight, especially at optical wavelengths. However, a 
current lack of instrumentation with resoved imaging po- 
larimety capabilities limits these possibilities. 
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